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ABSTRACT 

We have taken basis functions with cubic B-splines and weight functions with quartic B-splines in 
Petrov-Galerkin method to solve a boundary value problem of fourth order. In this method, the cubic B-splines and quartic 
B-splines are redefined into new sets of functions which contain the equal number of functions. To test the accuracy and 
efficiency of the method proposed, the numerical results obtained are presented in the form of absolute errors and found 
that the obtained results are giving a little absolute error. 
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INTRODUCTION 

Consider a general linear boundary value problem of fourth order 

Po 0)v (4) (0 + Px ( t)v"(t ) + p 2 (t)v'(t) + p 3 (0 v\t) + p 4 (t) v(t) = q(t) , a <t <b (1) 


Subject to the boundary conditions 


v(a) = A 0 , 

U° 

II 

§ 

, v\a) = A x , v\b) = C l 

(2a) 



or 


v(a) = A 0 , 

U° 

ii 

§ 

, v'\a) = A l , v"(b) = C x 

(2b) 



or 


v(a) = A 0 , 

v(b) = C 0 

, v\a) + <7 x v{a) =A, , v'(b) + <J 0 v(b) =C X 

(2c) 

where A 0 , C 0 , A b Ci, Oiand a 2 are real constants and p 0 ft), pi(t), p 2 (t), p 3 (t), p 4 (t) and q(t) 

functions defined in [a, b]. 

are continuous 


The fourth order boundary value problems arise in the areas of fluid mechanics, elasticity and quantum mechanics 
and in some allied areas of science and engineering. The existence and uniqueness for the solution of these problems are 
described in Agarwal [1]. The exact solutions of such type of boundary value problems are rarely available. Various 
numerical methods have been developed by many researchers [2-21]. Till now, the researchers did not use cubic B-splines 
as basis functions and quartic B-splines as weight functions in the Petrov-Galerkin method to solve fourth order boundary 
value problems of the type (l)-(2). That’s why we want to use the above method as our proposed method. 
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The subsequent sections are dealt with the justification of using Petrov-Galerkin method, a description of using 
proposed method with the types of boundary conditions (2), the procedure of solving the nodal parameters and the 
application of the proposed method on solving several examples of linear and nonlinear boundary value problems. 
By using quasilinearization technique [22], a nonlinear problem can be converted into a sequence of linear problems and 
the limit of solutions of these generated linear problems is the solution of the nonlinear problem. The conclusions are 
presented in the last section. 

JUSTIFICATION OF USING PROPOSED METHOD 

The approximate solution in Finite Element Method (FEM) can be obtained as a linear combination of basis 
functions which constitute a basis for the approximation space under consideration. Petrov-Galerkin method is one of the 
variational methods involved in FEM. The residual of approximation is made orthogonal to the weight functions in 
Petrov-Galerkin method. Regardless of properties of the differential operator defined in the given differential equation, a 
weak form of approximation solution for the differential equation exists and is unique under appropriate conditions [23, 
24], Further, if we pay sufficient attention to the boundary conditions [25], the weak solution tends to an exact solution of 
the differential equation. This means that the basis functions should become zero on the boundary where the essential 
(Dirichlet) type of boundary conditions are defined. Also the basis functions and the weight functions are equal in number. 

DESCRIPTION OF THE PROPOSED METHOD 

Cubic B -splines and Quartic B -splines 

The cubic B-splines and quartic B-splines are described in [26-28]. Space variable domain [a, b] is divided into 
spaced knots (which need not be spaced evenly) given by the partition a = fg < t\ < ... < 1 n _\ <t n = b . Six additional 
knots t. 3, t _ 2 , L i, t n+ 1, t n+ 2 and r n+3 are introduced which satisfy the relation 

f_3<f- 2 <f-i<fo and t n <t n+ i<t n+ 2 <t a+3 . 

Now the cubic B-splines Si(t)'s are defined by 


s,(t) 




x'(t) 


0 , 


U,-2> h+2 J 

otherwise 


Where 

i+2 

And 

r=i - 2 


Where {S.j(t), Soft), Sift),..., S n .j(t), S„(t), S„+i(t )} forms a basis for the space S 3 (^r) of cubic polynomial splines. 
Schoenberg [28] has shown that cubic B-splines are the unique nonzero splines of smallest compact support with the knots 
at t_3<f_2<t_i<to'^l < ^- • ■ 1 ^ tn+2< fi+3. 

In the same way, the quartic B-splines Ri(t)'s are defined by 
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[tj- 2 ’ ^,+ 3 ] 


otherwise 


i- 2 ’ l i + 3 


Where 



! + 3 


And 7 : ft) = ]^[ ft-t r ) 


Where {R- 2 (t), R.jft), Ro(t), Ri(t),..., R„.i(t), R n (t), R„ +1 (t)} forms a basis for the space S 4 (7r) of quartic polynomial 
splines by introducing two more additional knots t_ 4 ,t n+4 to the already existing knots t_ 3 to t n+3 . Schoenberg [28] has 


shown that quartic B-splines are the unique nonzero splines of smallest compact support with the knots at f_ 4 <f_ 3 <f_ 2 <f- 

I <t()<t| <■ . . 1 < t n <t t n+ <t M+ 2< f rl +3 7,1+4. • 


We define the approximation for v(t) as 




( 3 ) 


Where CC-S are the nodal parameters to be determined and Sj(t)'s are the cubic B-spline basis functions. In 

Petrov-Galerkin method, the basis functions should be zero on the boundary where the essential type of boundary 
conditions are prescribed. In the set of cubic B-splines {S.i(t), Soft), S](t),..., S n .:(t), S„(t), S n+ i(t)}, the basis functions S.ift), 
S 0 (t), S^t) do not become zero on the left boundary and S n _i(t), SJt) and S n+l (V do not become zero on the right boundary. 
So, it is necessary to redefine the basis functions into a new set of basis functions which become zero on the boundary. 

Applying the essential boundary conditions of (2), we get the approximate solution vf t) at the boundary points as 


= v fa) = v (t 0 ) = a_ l S_ l ft 0 ) + a Q S 0 ft 0 ) + a l S l ft 0 ) 


( 4 ) 



( 5 ) 


Eliminating CC x and OC n+x from the equations (3), (4) and (5), we get 


n 



(6) 


Where w ft ) = ^ S l ft ) + — ' — — S n+l ft) 
S-i(t 0 ) S n+l (t n ) 


( 7 ) 
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5,(f 0 ) 

5,(0 ' 

7 S_i(* 0 ) 

5,(0, 


5_ t (0, 7 = 0,1 

j - 2,3,...,n-2 


( 8 ) 


5 ) 

s j( 0- ' " ~ 5 n+1 (0, j = n-l,n 

s„M 


[Bjft), j — 0,1,..., «} are the new set of basis functions for the approximation v( t). Here wft) takes care of given set 

of essential boundary conditions and B ft) 's are vanishing at the boundary. In the proposed method, the new set of basis 

functions and weight functions should be equal in number. Here the number of basis functions in the approximation for vft) 
in (6) is n + 1 and the number of weight functions is n+4. So, it is necessary to redefine the weight functions into a new set 
of weight functions which are equal in number of the basis functions. 

Let us write the approximation for aft) as 


i =- 2 

Where R.ft)'s are the quartic B-splines. 

Method with Boundary Conditions (2a) 

Let us assume that approximation aft), given by (9), satisfies the conditions 
u(a) = 0,u(b) = 0,u' fa) = 0 

Using (9) and (10), we get the approximate solution for aft) at the boundary points as 
i 

ufa) = uft 0 ) = Yj o) = 0 

j=~ 2 


tfb)=uft n )= Y M(O = 0 


j=n-2 


i'fa) = uft Q )=Yj3 J R 'ft Q ) = 0 

j=~ 2 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


Eliminating , /? and B +1 from the equations (9) and (1 1) to (13), we get the approximation for aft ) ; 


uft) = (14) 

j = o 


where 




Tj(t) = < 

Pjft)~-y^—P-ift), 

J PM 

7 = 0,1 

(15) 


P M 

j = 2,3,4,...,n 
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Pjity- 


Rjfto) 

RM 


RXO- , k 2 ( 0 . 


*Vi+ 1 Vn ) 


j — 1 , 0,1 
7=2, 3,4,...,/? -3 
j — n — 2 , /? — 1 , n 
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(16) 


Now the new set of weight functions for the approximation u(t) is{73(f), j = 0 , 1 , ...,n} . Here 
Tj (t 0 ) = Tj (t n ) = T - (t 0 ) = 0 for all j . 

Applying the proposed method to (1) with the new set of basis functions { B ; (t), j = 0,1,...,/?} and with the new set 
of weight functions {73(f), j — 0,1,...,/?} defined in (15), we get 

t„ 

J toO)v <4) (0 + Pi (t)v"(t) + P 2 (t)v\t) + p 3 (t)v\t ) + P 4 (!)v(!)]T(!) <7i=j q(t)T(t)dt for i=0, 1 ,2, . . . ,n. (17) 


Integrating by parts the first two terms on the left hand side of (17) and after applying the boundary conditions 
mentioned in (2a), we get 


j Po {t)T t {t)v w {t)dt = Ci A 


dt 


dr 


(18) 


-} 4 t[a> (. 0 T i (t)]v\t)dt ~[p 0 C)T, (0], v"(f„ ) 


dt 


dt 


d 2 


j Pl (t)T {t)v"\t)dt = j (f ) 7} (t)]v\t)dt -^-[ P] (t)T (t)\ C, 


dt 


dt 


(19) 


Using (18), (19) and (6) in (17) and after rearrangement, we get a system of equations in the matrix form as 
Ka =f (20) 

where K = [&,-,] ; 


Kj = f { [- [Po m (0] ■ + [ Pl m (; t )] + p 3 (t)T (; t)]Bj ( t ) 
dt dr 

+p 2 (t)T i (t)B J "(t) + p 4 (t)T i (t)B j (t)}dt-^-[p 0 (t)T i (t)] Bj\t n ) 

dt 


for i=0, 1 ,2, . . . ,n; j=0, 1 ,2, . . . ,n. (21) 


/ = [/;]; 
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7 3 7 2 

f, = f [gw it) +[—j[ Po m (ty\- - -j [ a (0?; (0] - p 3 m ( t)]w\t ) 

■"o ru or 

-p 2 {t)T, (t)w"(t) - p 4 (t)T t (t) w(t) }dt - ^ [ p 0 (t)T f (()],„ c t 

+ ^ [Po(07;(0]l » A + f[ ;7o(07;(0 ]'.. VV " (? » )+ ^^ |( ^ (0 l C| 

and a = [a 0 a, ...aj r . 


Method with Boundary Conditions (2b) 

Let us assume that the approximation lift), given by (9), satisfies the conditions 
u(a)=0 ,u(b) = 0 ,u'\a) = 0 


( 22 ) 


(23) 


Using (9) and (23), we get the approximate solution for u(t) at the boundary points as 

u(a ) = u(t 0 ) = £ { 3 j R j (t 0 ) = 0 (24) 

;=- 2 


n + 1 

u(b) = u(t n )= X PjRj(t n ) = 0 (25) 

j=n-2 

u"(a) = u"(t 0 ) = ^ j3jR''(t 0 ) = 0 (26) 

7—2 

Eliminating /?_ 2 , /C and /?„+i from the equations (9) and (24) to (26), we get the approximation for u(t) as 

n 

u(t) = Y,PjTj(t) (27) 

7=0 


where 


P , (/„) 

/>(')- y=o,i 

(*<>) 

^•(0, j = 2,3,4,. 


(28) 


And {PXt),j = —1,0,1 ,..., 77 — l,/i} are defined in (16). 

Now {71(f), 7 = 0,1,..., n} are the new set of weight functions for the approximation u(t). Here 
Tj Q 0 ) = Tj (t n ) = T" (fo) = 0 for all j . 

Applying the proposed method to (1) with the new set of basis functions { B (t), 7 = 0,1,..., n} and with the new 
set of weight functions {71(f), j — 0 , defined in (28), we get 


Impact Factor (JCC): 2.6305 


NAAS Rating 3.19 



Numerical Solution of Fourth Order Boundary Value Problems by Petrov-Galerkin Method 
with Cubic B-splines as Basis Functions and Quartic B-splines as Weight Functions 


71 


<n t„ 

\\Po(t)v ( 4 } (t) + /?, (t)v"(t)+ P 2 (t)v"(t) + P 3 (t)v< \t) + P 4 (t)v(t)\ T(t) dt =jq(t)T(t)dt for i=0,l,2„ . . ,n. (29) 

Integrating by parts the first two terms on the left hand side of (29) and after applying the boundary conditions 
prescribed in (2b), we get 


—{PoWiit)) 

dt 


C,+ 


J P 0 {t)v w {t)T i {t)dt = - 

h 

t n 

} Pi(t)T i (t) V '"(t)dt = -j ^-[p^mi'xodt 


-r(p a d)m) 

dt 


d 2 


A ] + \—j(p 0 (t)T i (t))v"(l)dl 
J at 


r o h 


(30) 


(31) 


Using (30), (31) and (6) in (29) and after rearrangement, we get a system of equations in the matrix form as 
Ka =/ (32) 

where K = [k,j ] ; 


^ = r^TT^o(07;(0]-^[A(07;(0]+^(07:(0]s/ , (o 

rt nt 


dt 


dt 


for i=0, 1 ,2, . . . ,n; j=0, 1 ,2, . . . ,n. 


(33) 


+P 3 (t)T i (t)B j (f) + p 4 (f)7;.(f)5;(f)}Jf 


f = m 

1 2 i 

/;. = f t "{^(07;(o+[- T y[Po(07:(0]+-[A(07:(0]- J p 2 (07:(0]w"(o 
J, o dt dt 

-p 3 (07:(ovv / (o- J P4(07:(Ow(o}^+4[Po(07:(0] f wwia 

dt " dt 


(34) 


for i=0, 1, 2, ..., n 

and a = [a Q a x ...a n f . 

Method with Boundary Conditions (2c) 

Let us assume that the approximation u(t), given by (9), satisfies the conditions 

u(a)=0,u(b) = 0,u'(a) = 0 (35) 

Let us proceed as in the case of method with boundary conditions (2a). {T.(t), j — 0, 1 , tl) , as defined in (15), 
are the new set of weight functions for the approximation u(t). 

Applying the Petrov-Galerkin method to (1) with the new set of basis functions { Bj(t ), j = 0,1 n} and with 

the new set of weight functions {73 ( t ), j = 0,1,..., n} defined in (15), we get 
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}[p 0 (Ov <4) (0+ Pi(t)v"(t)+p 2 (t)v"(t)+p 3 (t)v< (t) +p 4 OM)fc(0 dt=jq(l)T(l)dl 


for i=0,l,2,...,n. (36) 


Integrating by parts the first two terms on the left hand side of (36) and after applying the boundary conditions 
prescribed in (2c), we get 


j p 0 W(t)v™(.t)dt = (C, -< 5l C 0 )-±-[p 0 {t)T l {t)\ (A -a,A) 


dt 


dr 


(37) 


’n J3 


J -^p 0 (t)T i (t)]v'(t)dt-—[p 0 (t)T i (t)]^ v"(t n ) 


dt 


l 1 J 2 


j Pi(t)T i (t)v"\t)dt^ j^[p l (t)T i (t)]v'(t)dt^[ Pl (t)T i (t)]^ (C,-a 2 C 0 ) 


(38) 


Using (37), (38) and (6) in (36) and after rearrangement, we get a system of equations in the matrix form as 
Ka=f (39) 

where K = [£,,] ; 


^ = f f '{[--j T [/ ? o(07;(0]+-7T[AW(0]+P3W(0]5/(f) 

Jt o dt dr 

+p 2 (07;.(05/ , (0 + P4(07:(05y(0}^-^-[Po(07;(OL Bj\t n ) 

dt 

for i=0,l,2,..., n; j=0,l,2,..., n. 

fi = f {9(071(0 +[^[p 0 (t)T i (t)\-^[p l (t)T i (t)\ - ^(OTKOJw'a) - P 2 (071(0 w"(0 

(if dt 

-p 4 (Or,(OH<O)*-|T[;'„(»)7:(O],,(c 1 -<7,c 0 )+|y[p 0 (O7’«)], t (A-^A) 
+|r[Po«)f(()], (c, -<7 2 c 0 )- 

for i=0, 1, ..., n 

and <2 = [a o a, ...Oj 7 . 

PROCEDURE OF SOLVING THE NODAL PARAMETERS 

n — 1 

A general element in the matrix A is given by Z'. , where I m = J + ‘ H.(f)r.(f)M(f)(if , r (0 are the 

m = 0 


(40) 


(41) 


cubic 


B-spline basis functions or their derivatives andw ; (0 are the quartic B-spline weight functions or their derivatives. Here 
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I m = 0 'f (f ,_2 > f ,+3 ) if j - 2 > tj +2 ) Cl (t m > ) = 0 . For the evaluation of each I m , we have used 4-point Gauss- 

Legendre quadrature formula. Due to this, the stiffness matrix A - is a eight diagonal band matrix. Solving the system Ka = / 
by using the band matrix solution package, we get the nodal parameter vector a. We have used the FORTRAN-90 code to 
solve the boundary value problems (1) - (2) by the proposed method. 

NUMERICAL EXAMPLES 

To test the accuracy and efficiency of the developed method, we solved five linear and four nonlinear fourth order 
boundary value problems. The obtained numerical results for each problem are presented in tabular forms. 

Example 1 

Consider the following linear boundary value problem 

V (4) + 4v = 1, - 1 < t < 1 (42) 


Subject to v(— 1) = v(l) = 0 , v\— 1) = 


sinh2-sin2 
4(cosh2 + cos 2) ’ 


sin2-sinh2 

v (1) = . 

4(cosh 2 + cos 2) 


The exact solution for the above problem is 


v = .25 


1-2 


sinh 1 sin 1 sinh t sin t + cosh 1 cos 1 cosh t cos t 
(cos 2 + cosh 2) 


Dividing the domain [-1, 1] into 10 equal subintervals, the numerical results obtained for this problem are 


presented in Table 1. 


Table 1: Numerical Results for Example 1 


t 

Absolute Error by 
Proposed Method 

-0.8 

8.095056E-06 

-0.6 

5.215406E-06 

-0.4 

7.696450E-06 

-0.2 

5.871058E-06 

0.0 

7.390976E-06 

0.2 

5.863607E-06 

0.4 

7.711351E-06 

0.6 

5.27501 IE-06 

0.8 

7.953495E-06 


Example 2 

Consider the following linear boundary value problem 

V (4) +tV = -(S + lt + t 3 )e‘, 0<t<l (43) 

Subject to v(0) = v(l) = 0, v'(0) = 1, v'(\ ) = —e. 

The exact solution for the above problem is V = f(l — t )e' . 


www.iaset.us 


editor @iaset.us 


74 


K. N. S. Kasi Viswanadham & S V Kiraninavi.Ch 


Dividing the domain [0, 1] into 10 equal subintervals, the numerical results obtained for this problem are 
presented in Table 2. 

Example 3 

Consider the following linear boundary value problem 

-v"-v = e t (t- 3), 0<t<l (44) 


Subject to v(0) = 1, v(l) = 0, V^O) = 0, /(l) = —e. 


The exact solution for the above problem is v = e (1 — t) . 


Dividing the 


domain [0, 1] into 10 equal subintervals, the numerical results obtained for this problem are 


presented in Table 3. 


Table 2: Numerical Results for Example 2 


t 

Absolute Error by 
Proposed Method 

0.1 

1.7665 3 3E-05 

0.2 

1.606345E-05 

0.3 

2.044439E-05 

0.4 

1.522899E-05 

0.5 

2.157688E-05 

0.6 

1.376867E-05 

0.7 

2.276897E-05 

0.8 

1.284480E-05 

0.9 

2.232194E-05 


Table 3: Numerical Results for Example 3 


t 

Absolute Error by 
Proposed Method 

0.1 

3.695488E-06 

0.2 

5.424023E-06 

0.3 

3.099442E-06 

0.4 

6.198883E-06 

0.5 

3.218651E-06 

0.6 

4.947 186E-06 

0.7 

4.649 162E-06 

0.8 

3.665686E-06 

0.9 

5.483627E-06 


Example 4 


Consider the following linear boundary value problem 


-v = — 4(2f cos? + 3 sin?), 0<f<l (45) 

Subject to v(0) = 0 , v(l) = 0 , v"(0) = 0 , v"(l) = 2 sin 1 + 4 cos 1 . 


The exact solution for the above problem is V = 


0 t 2 


- 1) sin t . 
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Dividing the domain [0, 1] into 10 equal subintervals, the numerical results obtained for this problem are 
presented in Table 4. 

Example 5 

Consider the following linear boundary value problem 

v (4) -3601v" + 3600v = -l + 1800f 2 , 0<f<l (46) 

Subject to v(0) = 1, v(l) = 1.5 + sinh(l), v'(0) - v(0) = 0, v'(l) - v(l) = -0.5 + cosh(l) - sinh(l) . 


The exact solution for the above problem is V = 1 + 0.5 t + sinh(f) . 


Dividing the 


domain [0, 1] into 10 equal subintervals, the numerical results obtained for this 


problem are 


presented in Table 5. 


Table 4: Numerical Results for Example 4 


t 

Absolute Error by 
Proposed Method 

0.1 

2.6971 10E-06 

0.2 

4.470348E-07 

0.3 

4.142523E-06 

0.4 

3.278255E-07 

0.5 

4.440546E-06 

0.6 

3.278255E-07 

0.7 

4.50015 IE-06 

0.8 

3.576279E-07 

0.9 

4.261732E-06 


Table 5: Numerical Results for Example 5 


t 

Absolute Error by 
Proposed Method 

0.1 

2.384186E-06 

0.2 

2.264977E-06 

0.3 

4.768372E-06 

0.4 

2.741814E-06 

0.5 

2.980232E-06 

0.6 

4.768372E-07 

0.7 

7.152557E-07 

0.8 

1.668930E-06 

0.9 

2.384186E-07 


Example 6 

Consider the following nonlinear boundary value problem 

v (4) -6e~ 4v = -12(1 + 0” 4 , 0 < f < 1 (47) 

Subject to v( 0) = 0, v(l ) = In 2, v'(0) =1, v'(l) = 0.5. 


The exact solution for the above problem is v = In(l+t). 


By using quasilinearization technique [22], the nonlinear boundary value problem (47) is converted into a 
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sequence of linear boundary value problems as 

v (nii) + [ 24 e _4v< " ) ]V(„ + i) = -12(1 + 0 ~ 4 +e ~ 4v< " ) [6 + 24 v ( b) ], n = 0,1, 2,3,... (48) 

Subject to (0) — 0, (1) — In 2, j (0) —1, v^+i) (1) — 0.5 . 


Here v (n+ i, is the (n+l) ,h approximation for v. The domain [0, 1] is divided into 10 equal subintervals and the 
proposed method is applied to the sequence of linear problems (48). The obtained numerical results for this problem are 
presented in Table 6. 


Table 6: Numerical Results for Example 6 


t 

Absolute Error by 
proposed Method 

0.1 

6.929040E-07 

0.2 

2.384186E-07 

0.3 

5.066395E-07 

0.4 

0.000000E+00 

0.5 

3.874302E-07 

0.6 

2.384186E-07 

0.7 

3.576279E-07 

0.8 

1.788139E-07 

0.9 

0.000000E+00 


Example 7 

Consider the following nonlinear boundary value problem 

v (4) - v /2 + vv" = -4r +e'(l + t 2 -4t), 0<t<l (49) 


Subject to v(0) = 1, v(l) = 1 + e, v'(0) = 1, v'(l) = 2 + e . 


2 t 

The exact solution for the above problem is V = t ” + e . 

By using quasilinearization technique [22], the nonlinear boundary value problem (49) is converted into a 
sequence of linear boundary value problems as 

2 

V(„|i)+V( n )V(" +1 )-2v(„)V( n+ i)+V(")V(„ +1 ) =-4C'+e f (l + f 2 ~40-(v( n )) +v (n) V(") , n=0,l,2„... (50) 

Subject to v (B+1) (0) = 1, v (B+1) (1) = 1 + e, v' B+1) (0) = 1, v' B+1) (1) = 2 + e . 


Here is the (tl + \) ,h approximation for v. The domain [0, 1] is divided into 10 equal subintervals and the 

proposed method is applied to the sequence of linear problems (50). The obtained numerical results for this problem are 
presented in Table 7. 
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Table 7: Numerical Results for Example 7 


t 

Absolute Error by 
Proposed Method 

0.1 

9.536743E-07 

0.2 

2.980232E-06 

0.3 

2.741814E-06 

0.4 

5.602837E-06 

0.5 

3.933907E-06 

0.6 

6.675720E-06 

0.7 

2.145767E-06 

0.8 

2.384186E-06 

0.9 

2.145767E-06 


Example 8 


Consider the following nonlinear boundary value problem 

= sinf + sin 2 l-fy"] 2 , 0<7<1 (51) 

Subject to v(0) = 0, v(l) = sinl,v"(0) = 0, v"(Y) = - sin 1 . 


The exact solution for the above problem is V = sin t . 

By using quasilinearization technique [22], the nonlinear boundary value problem (51) is converted into a 
sequence of linear boundary value problems as 


V(„ 4 !i) + [2 v") ] v ( : +1) = sin t + sin 2 1 + [v' n) f , n = 0, 1, 2, 3, ... 


(52) 


Subject to v (B+I) (0) = 0, v (B+1) (1) = sin 1 , v" +1) (0) = 0, v" +1) (1) = - sin 1 . 


Here is the (n +l) ,,, approximation for v. The domain [0, 1] is divided into 10 equal subintervals and 

the proposed method is applied to the sequence of linear problems (52). Numerical results for this problem are presented in 
Table 8. 


Table 8: Numerical Results for Example 8 


t 

Absolute Error by 
Proposed Method 

0.1 

8.456409E-06 

0.2 

1.586974E-05 

0.3 

2.205372E-05 

0.4 

2.57492 IE-05 

0.5 

2.6971 10E-05 

0.6 

2.521276E-05 

0.7 

2.157688E-05 

0.8 

1.531839E-05 

0.9 

7.987022E-06 


Example 9 

Consider the following nonlinear boundary value problem 
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v (4) -v 2 = e‘ —t 4 - e 2t - 2f V , 0<f<l 

Subject to v(0) = 1, v(l) = 1 + e, v'(0) - v(0) = 0, v'(\) - v(l) = 1 . 


(53) 


2 , J 


The exact solution for the above problem is V = t + e . 

By using quasilinearization technique [22], the nonlinear boundary value problem (53) is converted into a 
sequence of linear boundary value problems as 


v, <4 '-2v m v, =e' -r 4 -e 2 ' -2t 2 e' -v 2 M , n=0, 1,2, 3,... 

(«+l) («) (n+1) (n) ’ ’ ’ 

Subject to v (n+V) (0) = 1, v (B+1) (1) = 1 + v' B+1) (0) - v (B+1) (0) = 0, v' B+1) (1) - v (B+1) (1) = 1 . 


(54) 


Here V( n +]j is the (n + Y) ,h approximation for v. The domain [0, 1] is divided into 10 equal subintervals and the 

proposed method is applied to the sequence of linear problems (54). The obtained numerical results for this problem are 
presented in Table 9. 


Table 9: Numerical results for Example 9 


t 

Absolute Error by 
Proposed Method 

0.1 

1.668930E-06 

0.2 

1.192093E-07 

0.3 

2.026558E-06 

0.4 

1.668930E-06 

0.5 

3.695488E-06 

0.6 

2.145767E-06 

0.7 

5.245209E-06 

0.8 

3.814697E-06 

0.9 

4.0531 16E-06 


CONCLUSIONS 

In this paper, we have solved a general fourth order two point boundary value problem with three different cases 
of boundary conditions by the proposed method with cubic B-splines as basis functions and quartic B-splines as weight 
functions. The cubic B-splines and quartic B-splines are redefined into new sets of functions which contain the equal 
number of functions. To test the accuracy and efficiency of the developed method, it has been tested on five linear and four 
nonlinear fourth order boundary value problems. It is found that the obtained results are giving a little error. The strength 
of the developed method lies in the easiness of its application, accuracy and efficiency. 
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